arXiv:1504.04952v2 [cond-mat.mes-hall] 27 Dec 2015 


Tuning eqnilibration of qnantum Hall edge states in graphene role of crossed electric 

and magnetic fields 

Sudipta Dubej0 and Mandar M. Deshmuktj^ 

Department of Condensed Matter Physics and Materials Science, 

Tata Institute of Fundamental Research, Homi Bhabha Road, Mumbai 400005, India 

(Dated: December 29, 2015) 

We probe quantum Hall effect in a tnnable 1-D lateral snperlattice (SL) in graphene created using 
electrostatic gates. Lack of equilibration is observed along edge states formed by electrostatic gates 
inside the snperlattice. We create strong local electric field at the interface of regions of different 
charge densities. Crossed electric and magnetic fields modify the wavefunction of the Landau Levels 
(LLs) - a phenomenon unique to graphene. In the region of copropagating electrons and holes at 
the interface, the electric field is high enough to modify the Landau levels resulting in increased 
scattering that tunes equilibration of edge states and this results in large longitudinal resistance. 


Magnetotransport across one-dimensional snperlattice 
(SL) had been studied in two-dimensional electron gas in 
semiconductor heterostructures (2DEGS) [IHS], report¬ 
ing dissipationless transport across high potential bar¬ 
riers [T] and magnetic commensurability oscillations in 
longitudinal resistance [3]. The motivation was to study 
various competing length scales and energy scales be¬ 
tween tunable SL potential and quantum Hall system. 
Graphene offers the advantage of large cyclotron gap al¬ 
lowing quantum Hall effect to be observed at room tem¬ 
perature [THS]. Substrate induced SL in graphene in the 
presence of magnetic field led to the experimental ob¬ 
servation of Hofstadter butterfly physics [lUHIS- The 
ability to create abrupt 10 nm) tunable barriers in 
graphene allows new aspects to be explored. In addition, 
new physics, due to the role of crossed electric and mag¬ 
netic field, that cannot be seen in conventional 2DEGS 
can be studied in SL structures based on graphene. 

In this letter, we study magneto transport in an elec¬ 
trostatically defined ID lateral SL in graphene [T^. In 
our device we apply a perpendicular magnetic field and 
periodically modulate the charge carrier density in ad¬ 
jacent “ribbons” of graphene, tuning from an array of 
p-p’ (or n-n’) to an array of p-n’ junctions. Ghanging 
the magnetic field allows us to vary Ib relative to A; and 
changing the gate voltage allows us to tune the SL poten¬ 
tial strength relative to LL spacing. The relative abrupt¬ 
ness, bipolarity of charge carriers, large modulation and 
unequally spaced LLs distinguishes the present work from 
the previous work on ID SL using 2DEGS systems P- 

HHl]. 

Apart from the length scales, we also study the energy 
scales involved. The competition between SL amplitude 
(Vb) and LL spacing (Jiujc, where h = h/2'!T, h being 
the Planck’s constant, and Wc is the cyclotron frequency) 
gives rise to three regimes. When Vq » hujc, SL effect 
dominates giving rise to extra Dirac points m- In the 
other extreme when Vq << huJcj quantum Hall effect in 
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graphene is restored |15j . However, the situation is more 
complex and little explored when Vg and fiUc have com¬ 
parable contribution, and we have experimentally probed 
this regime in graphene. 

The goal of our work is to extend quantum Hall stud¬ 
ies beyond single top-gate in graphene [TBHTS] . Our work 
is the first experimental report on magnetotransport in 
multiple top-gates on graphene and we probe the physics 
of equilibration along the narrow region in graphene de¬ 
fined electrostatically. Our main observation is that when 
Vq is comparatively small in the unipolar region, the edge 
states do not equilibrate along this narrow region de¬ 
fined electrostatically. The extent of equilibration can 
be tuned in the bipolar region where the electric field is 
relatively large. In this regime, electric field significantly 
modifies the Landau level wavefunctions, increasing scat¬ 
tering, which is reflected in increased equilibration and 
large longitudinal resistance. 

We create a ID tunable SL, of period A, by fabri¬ 
cating an array of thin finger gates on graphene. The 
schematic of a device is shown in Figure [^a), and Fig¬ 
ure shows false colored scanning electron microscope 
image (details of fabrication in Section I of Supplemental 
Material). The geometric width of each top-gates is ^ 
30 nm and they have a period of A = 150 nm. The effec¬ 
tive electrostatic width of the top-gates felt by the charge 
carriers in graphene is larger due to the finite thickness 
of the top-gate dielectric [T3] (details in Section VH of 
Supplemental Material). 

In our device, graphene consists of two alternating re¬ 
gions - one where the charge carrier density is controlled 
only by the back-gate (BG region); and the other where 
the charge carrier density is controlled by both the top- 
gate and the back-gate (TG region). The difference in 
charge carrier density between BG and TG regions gives 
rise to a SL whose amplitude (Vb) is controlled by Vbg and 

Vtg. (Vb = y/^hvF{sgn{CbgVbg)\J - sgn{CtgVtg + 

CbgVbg) ^ where Vtg {Vbg) is the top- 
gate (back-gate) voltage, Ctg {Cbg) is the capacitance 
per unit area of top-gate (back-gate), e is the electronic 
charge and vp is the Fermi velocity.) (Details of calcu- 
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FIG. 1. Device geometry and different regimes of edge state 
transport, (a) Schematic of a device, (b) False colored scan¬ 
ning electron microscope image of a device with a zoomed-in 
image of the finger like top gates, (c) Depiction of the peri¬ 
odic ID potential and the influence of the two gates in BG 
and TG regions, (d) Parameter space of Vtg and Vbg showing 
the different type of charge carriers in adjacent regions, (e) 
Line plot of resistance as function of Vi,g when Vtg is biased 
in charge neutral region at 0 T and 2 K. (f) Line plot of re¬ 
sistance at 0 T and 2 K when TG regions are n-doped. (g) 
Schematic depicting edge state transport in the unipolar re¬ 
gion with Vbg greater than vtg, and the bipolar region. The 
number of top-gates is four in this schematic. The BG region 
is denoted in purple and the TG region is shaded in blue. 


lation and plot of Vb as a function of Vbg and Vtg is in 
Section II of Supplemental Material.) 

We measure zero-bias four-probe longitudinal resis¬ 
tance (Rxx) while varying gate voltages at different mag¬ 
netic fields {B) at a temperature of 2 K. The charge neu¬ 
tral point is at Vtg = -0.1 V and Vbg = -2 V (Figure[^e)) 
suggesting low unintentional doping. The mean free path 
in our device is ~ 70 nm and phase coherence length is ~ 
600 nm at 2 K m- As B increases, the magnetic length 

{Ib = \J^s) decreases and the charge carriers encounter 
smaller periods of SL until they are confined within a sin¬ 
gle BG or TG region. Well resolved LLs start to appear 
only beyond 2 T (see Supplemental Material Section III 
for complete LL fan diagram). In this work we look at 
quantum Hall effect in ID SL at = 22 and vary Vq 
upto ^ 375 meV. 

The resistance depends on the filling factor v = nh/Be 
in the adjacent regions [M113 uni mi- At a constant 
magnetic field B, the filling factor in the adjacent regions 
is determined by Vbg and Vtg. Depending on Vbg and Vtg 
at a given B, we have same type of edge states in the 
unipolar region, and in the bipolar region, electron and 
hole edge states co-propagate along the junction (Fig- 
ure[^g)). In this device, we have 37 top-gates. So, if the 
edge states are formed under all the top-gates and they 
equilibrate along all the edges of the top-gates, then the 


four-probe resistance (Rxx) plateau is given by 


R 


XX 


h ^Iktgl-ktgll 

ktalkbgl 

-^(ktgl + Fbal) 

e2 \Vtg\\Vbg\ 


l^tgVbg > 0 (I) 

VtgVbg < 0 (II) 


where N is the number of top-gates which is 37 in our 
device. (Landau-Biittiker formalism to obtain Rxx for 
multiple top-gates in Section V of Supplemental Mate¬ 
rial.) 

Figure a) shows the colorscale plot of zero bias four- 
probe resistance as a function of Vtg and Vbg at 14 T. We 
observe diamond shaped regions in the parameter space 
that represent integer filling factors in adjacent regions 
set by Vtg and Vbg. The filling factors in the two alter¬ 
nating regions are indicated as {vtg,Vbg). Figure [^b, c) 
show line slices of Rxx as function of Vbg at vtg = 2, 6. 
(Line slices at vtg = 10 and 14 are presented in Section IV 
of Supplemental Material.) We do not observe large re¬ 
sistance as predicted by Equation I and 11. The Rxx is 
relatively high, and does not show a plateau, in the bipo¬ 
lar regime, but not as high as predicted by Equation 11. 

In the line slices (Eigurej^b, c)), the green curve is the 
experimental data. The black dashed lines correspond 
to the calculated plateau for V = 1 in Equation I and 
11. We find that in the unipolar regime, denoted by the 
blue region in Figure [^d), the potential Vq is small and 
the plateaus are well described by Rxx = 

6 l^tgW^bgl 

there is good agreement between measured experimental 
data and expected plateau values for a single top-gate in 
Figure [^b,c). 

Our experimental observation show plateaus corre¬ 
sponding to a single top-gate and not 37 top-gate, and so 
we do not have equilibration in our device inside the su¬ 
perlattice. The equilibration occurs only at the extreme 
edge of the top-gates near the voltage probes as illus¬ 
trated in the schematic in Figure ^d). In the schematic, 
the yellow probes are the real probes and the gray probes 
are the virtual probes denoting equilibration at that edge. 
The virtual voltage probes are used to calculate resis¬ 
tance using Landauer-Biittiker formalism |22j . One pos¬ 
sible reason for the lack of equilibration is that the edge 
states under the top-gate are defined electrostatically, 
where the potential varies smoothly due to the finite 
thickness of the top-gate dielectric [23]. Equilibration 
requires inter edge state scattering or ohmic contacts so 
that all the edge states are at the same chemical poten¬ 
tial, and this does not happen due to the short length 
along the physical edge of graphene [24H27|. 

In the bipolar region, plateau coincides with that ex¬ 
pected for a single top-gate when vtg = 2(-2) and Vbg = 
-2(2). When we have one edge state for both electrons 
and holes circulating in adjacent regions, the resistance 
plateau is seen at hje^ (25.8 kfl). The plateau at (vtg^Vbg) 
= (2,-2) results not only because the E is low in this state 
compared to any other state in the bipolar region but also 
because this is a special state where the charge carriers 
belong to the same LL (n = 0) and that the LL gap is 
maximum between n = 0 and n = 1 LL. Narrow TG and 
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FIG. 2. Longitudinal resistance (Rxx) at i? = 14 T when LLs 
are well resolved in three different case of edge state trans¬ 
port. (a) Rxx as function of Vtg and Vbg at 14 T at a tem¬ 
perature of 2 K. The white lines denote line of constant Vtg- 
Rxx as function of Vbg for (b) vtg = 2 and (c) Vtg = 6. The 
green curves in (b), (c) denote experimental value with the 
black dashed lines representing the resistance plateau in case 
of equilibration of edge states for a single top-gate. {{vtg,Vbg) 
values are denoted within the parenthesis.) (d) Schematic 
depicting equilibration along the extreme edges in unipolar 
region. The yellow probes are the real contacts and the gray 
probes are the virtual voltage probes. 
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FIG. 3. Effect of E on equilibration (a) Variation of LLs 
along the length of device when (vtg,Vbg) = (-2,6). The brown 
region in the schematic denotes the junction between TG and 
BG region where E exists, (b) Spatial variation of potential 
and E along the length of graphene in bipolar region (Vi,g = - 
30 V and Vtg = 3 V). Yellow regions denote geometric position 
of top-gates, (c) Effective magnetic length as function of E at 
B = 14 T. (d) Gontour of maximum E at 1.2x10^ V/m and 
1.4x 10^ V/m as function of Vi,g and Vtg overlaid on measured 
Rj,!, in bipolar regime at 14 T. 


BG regions inside the superlattice leads to leaky barrier 
and transmission of edge states resulting in a plateau 
corresponding to a single top-gate [T] . 

However, when neighboring regions have other edge 
states, that is, (vtg x v^g) < 0 and \vtg\ or \iyi,g\ is greater 
than 2, we find resistance significantly larger than h/e^ 
(for example, (2,-6) in Figure [^b), (6,-2) in Figure [^c)); 
a feature not generally seen in single p-n’-p junction in 
the quantum Hall state. The feature of plateaus in R^x 
corresponding to equilibration for a single top-gate in 
the unipolar region and deviation from this picture in 
the bipolar region, is quite robust and is also observed in 
another device with 35 top-gates (details in Section X of 
Supplemental Material). 

The large resistance with maximum Rxx 200 kH at 
14 T, is seen along the diagonal direction in Figure [^a); 
it is precisely in this diagonal direction of the parameter 
space that Vq increases. In the bipolar region, though 
the value of Rxx is higher than that for a single top- 
gate, it is not 37 times the value for a single top-gate, 
and thus there is no full equilibration. In this region the 
array of top-gates does not behave as a single top-gate. 
The resistance value depends on the gate voltages and 
the magnetic field, implying that the equilibration can 
be modihed by electric field. 

We now try to understand the reason for this obser¬ 
vation and note that Vq is relatively large in this re¬ 
gion. Figure]^ a) shows a LL diagram along the length 


of the device when {vtg,Vbg) = (—2,6); where we have 
co-propagating electron and hole edge states at the junc¬ 
tion. At the interface of these regions (shaded brown in 
Figure [^a)) there is an electric field (A), due to the SL, 
and it has a significant effect on the LL wavefunctions. 

To get an accurate idea of the magnitude of E in our 
devices we performed numerical simulation of the electro¬ 
statics using finite element method. At a given Vtg and 
V/g, the charge carrier density induced along the length 
of graphene is calculated, from which potential and E is 
obtained. (Details of calculation in Section VH of Sup¬ 
plemental Material.) Spatial variation of potential and 
E along length of graphene in bipolar region is shown in 
Figure ^b). At a given Vtg and V/g, the maximum E is 
obtains which is higher in the bipolar region compared 
to the unipolar region. 

We note that in the bipolar regime, Vq created is much 
larger than fioJc- Large Vq leads to large E in the re¬ 
gion between BG and TG region which modifies the LLs 
locally. It has been shown by Lukose et al. and 

later extended by Gu et al. [15] for the case of a top- 
gate geometry that the LL spectrum and the wavefunc¬ 
tions in crossed E and B are fundamentally modified in 
graphene, an aspect that is not observable in conven¬ 
tional 2DEGS semiconductors. LL wavefunction is mod¬ 
ified in two ways. Effective magnetic length in the pres¬ 
ence of E can be written as I'b = Ib/{1 - {E/vfB)Y^^ 
[28] . I'g increases with increasing E and rises rapidly 
when E approaches vpB as seen in Figure [^c). Sec- 
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ondly, the E mixes the Landau levels. 

Contours of maximum E at 1.2x10^ V/m and 
1.4x10^ V/m are overlaid on the measured data at 14 T 
as shown in Figurej^d). Figure|^d) shows that the con¬ 
tours he along the region where we have high resistance 
state and departure from the value for a single top-gate. 
We argue that in our device geometry, E created is high 
enough to modify LLs which is reflected in charge trans¬ 
port measurements. The spatial extent of the wavefunc- 
tion (I'g) increases with increasing E and approaches the 
width of TG (or BG) region. This leads to increased 
spatial overlap of the wavefunction within TG (or BG) 
region of the superlattice resulting in enhanced scatter¬ 
ing. Secondly, as function of E, there is LL mixing [28] 
which results in a non zero matrix element essential to 
cause scattering and equilibration. Let us now examine 
the various lengthscales of our system that support this 
scenario. From Figure [^b) we find that the effective 
electrostatic width of the top-gate for this configuration 
is ^50 nm and the extent of the region with high electric 
field (between 1.2 - 1.4x10^ V/m) is ^15 nm. In addi¬ 
tion, the magnetic length in the presence of transverse 
electric in this region is ~11 nm (see Figure [^c)). If one 
considers the scenario of (6,-2) state, we find that with 
three edge states with width 11 nm will have significant 
overlap with two adjacent regions and equilibration will 
be enhanced due to scattering. 

In the presence of high electric held, the resistance de¬ 
pends on E and B applied. For example, resistance val¬ 
ues at {vtgjVbg) = (6,-2) and (2,-6) at 14 T are different 
as the maximum E is different in the two regions (Plot of 
maximum E as a function of gate voltages in Section VII 
of Supplemental Material) — this strongly suggests elec¬ 
tric held tuning of equilibration. 

We think disorder does not play an important role. 


Because, in our sample, the disorder potential, estimated 
from the FWHM of the Dirac peak, is 71 meV [13], and 
thus is smaller than Vq in the bipolar region. So, the 
tuning of equilibration of edge states in our sample, which 
is seen at higher Vq, is due to the E at the interface of 
TG and BG regions. However, recent work of Kumada 
et al suggests that the disorder along the length of p- 
n junction could play an important role in equilibration 
[33 ED]- Further experimental and theoretical studies 
need to be carried out to probe the role of disorder in 
periodically modulated structures. 

Details of LL modihcation in unipolar region at lower 
B of 3.5 T in Section IX of Supplemental Material. 

Our experiments with tunable superlattices suggest 
that tuning of the equilibration of edge states in graphene 
can be done using E at interfaces, which cannot be real¬ 
ized in conventional 2DEGS. In addition, the nature of 
the state that emerges after the collapse of the LLs is 
little understood and possibility of existence of correla¬ 
tions has been speculated ESEI]. The close proximity of 
co-propagating electron and hole edge states can be used 
to construct large class of topological states [33] and also 
offers an opportunity to study excitonic effects, this has 
been recently explored in bilayer quantum Hall systems 
[33] . There have been predictions of correlated states in 
V = Q LL of graphene and LL mixing can enable explo¬ 
ration of such phases [53] . 
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lowship of Department of Science and Technology and 
Department of Atomic Energy of Government of India 
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I. DEVICE FABRICATION 

Graphene is mechanically exfoliated on degenerately doped silicon (Si^'*') substrate having 
300 nm of silicon dioxide. Si"*"^ acts as global back-gate, that is, it times the charge carrier 
density of the entire graphene flake. Standard e-beam lithography is used to pattern Cr/Au 
contacts in Hall bar geometry. The device is then thermally annealed to remove residual 
PMMA from graphene. 3 nm of AI2O3 is then deposited using e-beam lithography which 
acts as seed layer for subsequent deposition of 20 nm AI2O3 using atomic layer deposition. 
AI2O3 acts as top-gate dielectric above which narrow (~ 27 nm wide), periodic fingers of 
Pd are patterned using e-beam lithography. All the fingers of top-gate are connected to one 
electrode outside the graphene flake so that all the hngers are at the same potential. The 
period of these hngers is 150 nm. 


II. VARIATION OF AMPLITUDE OF SUPERLATTICE POTENTIAL 
WITH BACK-GATE AND TOP-GATE VOLTAGE 

In our device there are two alternating regions in graphene - one without a top-gate 

(BG region), where the charge carrier density is controlled only by the back-gate; and the 

other under a top-gate (TG region), where the charge carrier density is controlled by the 

back-gate and the top-gate. The amplitude of superlattice potential (Vq) created is given 

by the difference of doping in the BG and TG region. In case of no unintentional doping, 

the number density of charge carriers in the BG region is given as nhg= ^ ; and in the 

TG region, the net number density of charge carriers is the algebraic sum of the density of 

c v -\~c v 

carriers induced by the two gates ntg= ] where Cbg {Ctg) is the capacitance per 

unit area between graphene and back-gate (top-gate) and Ug (kjg) is the back-gate (top- 
gate) voltage. Using parallel plate capacitor geometry Cbg is given by e^erld, where cq = 
8.85x10^^ F/m, = 3.9 and d = 300 nm. Ctg is obtained from the slope of charge neutral 
line in the plot of resistance as function of Ug and Vtg at zero magnetic field (Figure SI (a)). 
From the slope, we get Ctg = ISxCbg- In case of graphene, due to its linear dispersion, 
Ep = ^/phvp^/n, where Vp is the Fermi velocity and Ep is the Fermi energy, with all 
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energies measured relative to the charge neutrality point^. Thus Vq is given by 

Fo = ^hvp (^gn{C,gV,g) - sgn{CtgVtg + (i) 

Magnitude of Vq is referred as Vq since it represents the amplitude of superlattice potential. 
Figure Sl(b) shows the variation of Vq with Vlg and Vtg. 



FIG. SI. Amplitude of superlattice potential (a) Resistance as function of Vtg and Vbg at zero 
magnetic field at temperature of 2 K. White arrow points to the charge neutral line, (b) Contourplot 
of Vo as a function of Vbg and Vtg. Vq is higher in the bipolar region where there is a series of p-n’ 
junctions compare to the unipolar region. 


III. QUANTUM HALL IN THE ABSENCE OF SUPERLATTICE 
POTENTIAL 

Figure S2(a) shows the magnitude of Hall resistance as function of Vbg and magnetic held 
(B) when Vtg is biased at charge neutral point {Vtg = 0.1). The fan diagram eorresponds 
to monolayer graphene. Figure S2(b) is a line plot at 10 T showing plateau at u = 2 {= 
12.9 kH), u = 6 {= 4.3 kH) and u = IQ {= 2.5 kH), where u is Filing factor. We can observe 
the Landau levels from 2 T. 

IV. OBSERVED RESISTANCE FOR DIFFERENT FILLING FACTORS IN 
TG AND BG REGION 

The Filing factor is given by u = nh/Be, where n is the charge carrier density, h is the 
Planck’s constant, B is the magnetic Feld and e is the electronic charge. For constant R, 
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FIG. S2. (a) Magnitude of Hall resistance as a function of Vbg and B (b) Line plot of magnitude 
of Hall resistance as function of Vbg at 10 T showing plateaus at i/ = ± 2, ± 6, + 10. 


filling factor in BG region {ubg) is controlled by Vbg, and filling factor in TG region (utg) is 
controlled by both Vtg and Vbg. 

Line slices showing longitudinal resistance (Rxx) as function of Ubg for Utg = 2 and 6 are 
shown in Figure 2 in main text. For Utg = 10 and 14, we measured Rxx as function of Vtg 
and Vbg (Figure S3(a)) at a lower held of 8 T. Figure S3(b,c) shows Rxx as function of Ubg at 
Utg = 10 and 14. In the line slices (Figure S3(b, c)), the green curve is the experimental data 
with the black dashed line marking the calculated plateau in case of edge state equilibration 
for a single top-gate. Plateaus coincide with the black dashed lines in the unipolar region 
but deviates from them in the bipolar region. 
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FIG. S3, (a) Longitudinal resistance (Rxx) as function of Vtg and Vbg at 8 T. Rxx as function of 
Ubg for (b) Utg = 10 and (c) utg = 14. The green curves in (b), (c) denote experimental value with 
the black dashed lines representing the resistance plateau in case of equilibration of edge states for 
a single top-gate. {{utg,Ubg) values are denoted within the parenthesis.) 
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V. RESISTANCE IN CASE OF EQUILIBRATION OF EDGE STATES 

The schematic of our device geometry is shown in Figure S4. The yellow probes denote 
real contacts used to measure the resistances. The TG regions are denoted in blue and BG 
regions in purple. 



FIG. S4. Schematic showing our device geometry with multiple top-gates. In this schematic, six 
top-gates are shown. Yellow probes denote real contacts. Longitudinal resistance is measured 
between voltage terminal 2 and 3, and Hall resistance between voltage terminal 2 and 6. 


We now describe the resistance observed when ballistic conduction occurs via edge states 
using the Landau-Biittiker formalism^. We present the calculation for two top-gates since the 
matrix size increases with increasing number of top-gates. The procedure can be generalized 
for N top-gates, where N is 37 in our device geometry. The schematic (Figure S5)(a) 
shows two top-gates and we assume equilibration among edge states. The assumption of 
equilibration is incorporated in our calculation using virtual voltage probe which enforces 
all the edge states to be at the same chemical potential in that terminal. The real voltage 
probes are marked in yellow. The gray probes denote virtual voltage terminal which enforces 
equilibration along the electrostatic edge of the top-gates. \ui,g\ is greater than \utg\ and the 
adjacent region have same polarity of charge carriers. By counting the number of current 
carrying channels that start from terminal p and end in terminal q, transmission function 
Tpg is obtained. The conductance matrix Gpg is given by multiplying the following matrix 
with e^/h. 
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FIG. S5. Schematic showing two top-gates with the yellow probes as the actual contacts used to 
measure Rxx and Rxy- Gray probes are the virtual voltage terminals used to enforce equilibration 
of the edge states along the electrostatic edge of the top-gate, (a) Equilibration when \i'hg\ is 
greater than {utgl in the unipolar region, (b) Equilibration along all the electrostatic edges of the 
two top-gates, with \i'hg\ and {utg] having opposite polarity in the bipolar region. 
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1 2 3 4 5 6 7 8 9 10 11 12 13 14 

100 0 0 0 0000 0 0 0 0 \ubg\ 

2 \ubg\ 0 0 0 0 0000 0 0 0 00 

3 0 \ng\ 0 0 0 0000 0 o o oo 

4 0 0 \i^tg\ 0 0 0000 0 0 \ubg\-\utg\ 0 0 

5 0 0 0 [Ubgl 0 0000 0 0 0 00 

6 0 0 0 0 \vtg\ 0 0 0 0 {I'bgl — Wtg\ 0 0 0 0 

7 0 0 0 0 0 \ubg\ 000 0 0 0 00 

8 0 0 0 0 0 0 [ubgl 0 0 0 0 0 0 0 

900 0 0 0 00 \ubg\ 0 0 0 0 0 0 

10 00 0 0 0 000 \ubg\ 0 0 0 0 0 

11 0 0 0 0 \iyhg\ — l^tgl 0 0 0 0 \utg\ 0 0 0 0 

12 00 0 0 0 0000 0 \ubg\ 0 0 0 

13 0 0 \Ubg\ — \Utg\ 0 0 0000 0 0 \utg\ 0 0 

14 0 0 0 0 0 0 0 0 0 0 0 0 l^sl 0 

( 2 ) 

The current at a terminal p is given by Ip = Gpq{Vp — Vg). All of these equations are 
not independent, so, voltage at the drain terminal (terminal 8) is chosen to be zero, and we 
omit the row and the column corresponding to this terminal. Inverting the resulting matrix, 
voltages at the terminals are obtained. Also, the current at the voltage terminals are zero. 

The two probe voltage {V 2 -probe) is given by (lA — hg)) the four-probe longitudinal voltage 
{Vxx) is given by (lA — Vy) and the Hall voltage (Vxy) is given by (lA — Vu)- Dividing by 
current Ii, the resistances are obtained. Using the above formalism. 





Similarly, for \utg\ greater than \uhg\ for N top-gates, 


N\i^tg\ - (N - l)\ubg\ 


Rx 


Wtg\- Wbg\ 

h \T^bg\\ytg\ 




xy 


1 

h \ybgy 


( 5 ) 


The schematic for electron edge states and hole edge states in the adjacent region, that 
is, \vtg\ and \ubg\ having opposite polarity for two top-gates is shown in Fignre S5(b). The 
virtual gray terminal ensures equilibration among electron and hole edge states. The above 
mentioned procedure can be followed to obtain the resistance as, 


, _ {N + ^Wbg\ 

WKi 


Rx 


Wtg\ + Wbg\ 

h Hg\Wtg\ 


R 


■xy 


1 

h \Vbgy 


( 6 ) 


where N is the number of top-gates. 

In the unipolar region, the observed Rxx for 37 top-gates corresponds to the equilibration 
in case of a single top-gate. One possible reason is the narrow region in graphene dehned 
electrostatically which facilitates the transmission of edge states along the length of graphene 
without equilibration with locally circulating edges. The schematic of this case is shown in 
Figure S6 where the gray probes along the extreme edges denote equilibration along those 
two particular edge. Another reason could be that the extreme edges of the top-gate are 
near (~ 300 nm) the actual voltage probe which helps in the equilibration along those edges. 



FIG. S6. Schematic with yellow probes as the real contacts and the gray probes as the virtual 
voltage terminals denoting equilibration of edge states along that edge. Schematic showing equili¬ 
bration along the extreme edges of the top-gate in the unipolar region where {I'bgl is greater than 
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VI. VARIATION OF HALL RESISTANCE WITH GATE VOLTAGES 

Lack of equilibration or deviation from Landau-Biittiker picture is also observed in the 
Hall resistance [Rxy]- In our device geometry, the electrodes measuring R^y is outside the 
array of huger gates and hence the edge state equilibration gives R^y to be a function of 
only Vbg as seen in the previous section. 

Figure S7(a) shows Rxy measured as function of Vbg and Vtg at 14 T. A horizontal line slice 
showing plateaus at Vbg = ±2, ±6 is shown in Figure S7(b). Figure S7(c,d) shows Rxy as 
function of Vtg for Vbg = 2, 6. Rxy being a function of only Vbg should be independent of Vtg. 
We observe that when there are same type of charge carrier in the adjacent region, that is, 
for positive Ubg (positive Vtg), Rxy is almost constant. However, when we have electrons and 
holes in the adjacent region (positive Ubg and negative Vtg), Rxy does not remain constant. 
As we go from p-p’ (or n-n’) to p-n’, Rxy, from its constant value shows a dip and then 
increases toward the constant value. This further supports our picture that there is enhanced 
scattering in the elements of the junction. 


VH. NUMERICAL SIMULATION OF ELECTROSTATICS USING 
COMSOL 

We have modeled our device geometry in COMSOL software and obtained the electric 
field from numerical simulation of the electrostatics. In the simulation, we used 15 top-gates 
along the length of graphene and obtained the induced charge carrier density in graphene 
for a given back-gate and top-gate voltage. This induced charge carrier density is periodic 
along the length of graphene from which the potential is obtained. The electric held in 
the plane of graphene is obtained by taking the numerical derivative of the potential. The 
potential and the electric held variation along length of graphene for Vbg = -30 V and Vtg 
= 3 V in the bipolar region are shown in Figure S8(a,b). In the simulation, the geometrical 
width of the top-gate electrode is 30 nm, and the top-gates are 150 nm apart, similar to 
our device geometry. Due to the hnite thickness of top-gate dielectric, there is a spread of 
the electric held lines from the top-gates. So, electrostatically, the ehective width of the 
top-gated region in graphene is larger than 30 nm. For clarity, magnihed images with 3 
top-gates along length of graphene are shown in unipolar region (Figure S8(c)) and bipolar 
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FIG. S7. Transverse resistance (Rxy) at 14 T at a temperature of 2 K. (a) Magnitude of Rxy as a 
function of Vtg and Vbg. (b) Magnitude of Rxy as a function of Vbg when Vtg =-0.1 V. Plateaus at 
Vhg = ±2 and ±6 are observed. Rxy as function of Vtg at (c) Vbg = 2 and {d)i'bg = 6. In our device 
geometry, electrodes measuring Rxy is outside the array of finger gates and hence Rxy should be 
a function of only i^bg, that is, it should be constant with Vtg. However, Rxy decreases from its 
constant value when Vtg is such that there is opposite type of charge carrier in adjacent region. 


region (Figure S8(d)). We observe the magnitude of electric field to be larger in bipolar 
region compared to the unipolar region. Another interesting thing to note is that not only 
the magnitude but also the shape of the potential is different in the two cases. 

Numerical simulation of spatial variation of electric held along length and width of 
graphene in bipolar region {Vbg = -30 V and Vtg = 3 V) is shown in Figure S9(a). Electric 
held varies periodically along the length of graphene with maxima and minima at the edge 
of the top-gate, and is constant along the width of graphene. From the simulation, the max- 
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FIG. S8. Numerical simulation with 15 top-gates along length of graphene using COMSOL. (a) 
Variation of potential along length of graphene when Vtg = 3 V and Vbg = -30 V. (b) Spatial 
variation of electric field when Vtg = 3 V and V^g = -30 V. Magnified image of spatial variation of 
potential and electric field in (c) unipolar region when Vtg = 3 V and Vbg = 30 V and (d) bipolar 
region when Vtg = 3 V and Vbg = -30 V. Yellow regions denote geometric position of top-gates. 

imum electric field for a given Vbg and Vtg is obtained. Figure S9(b) shows the maximum 
value of electric field as a function of Vbg and Vtg in the bipolar region. 

VIII. EFFECT OF ELECTRIC FIELD ON EQUILIBRATION 

As seen in the previous section, the potential and the electric field is high in bipolar 
region where we have a series of p-n’ junctions compared to unipolar regime where there 
is series of n-n’ (or p-p’) junctions. In the unipolar region, there is lack of equilibration 
of edge states along the electrostatic edge and the array of top-gates behave as a single 
top-gate. The resistance plateaus depends only on the combination of filling factors, that is, 
resistance plateau at {i^tg^^bg) = (2,6) and (6,2) are same. In the bipolar region, the electric 
field modifies this lack of equilibration along the electrostatic edge such that the resistance 
plateau does not depend only on the combination of Filing factor in the adjacent region. The 
electric field causes modification of Landau level wavefunction. With increasing electric held 
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Length along graphene (nm) ^bg 


FIG. S9. Numerical simulation of electric field in our device in the bipolar region, (a) Spatial 
variation of the electric field along the length and width of graphene when Vtg = 3 V and Vhg = 
-30 V. The effective electrostatic widths of the TG and the BG regions for this gate voltage are 
marked in white, (b) Maximum electric held as function of Vtg and Vhg- 


the spread of the wavefunction increases such that it is of the order of the width of TG or BG 
region, resulting in enhanced scattering across Landau levels. This effect is manifested in 
charge transport when electric field becomes comparable to magnetic field times the Fermi 
velocity. Variation of effective magnetic length at 8 T as a function of gate voltages is shown 
in Figure SIO. 



FIG. SIO. Effective magnetic length at 8 T as a function of Vtg and Vbg. The effective magnetic 
length is calculated from the maximum electric held obtained from the simulation as shown in the 
previous section. 


In the bipolar region, the value of resistance depends on the electric and magnetic held 
in that region. For example, resistance value at {utg,Ubg) = (2,-6) and (6,-2) are different as 
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the potential and the electric field are different in these two regions. Similarly resistance 
value for {vtg,i^bg) = (6,-2) at 14 T and 8 T are different as shown in Figure Sll(a,b). 



Vbg 



FIG. Sll. Dependence of magnetotransport on electric field. Rxx as function of Vf^g for vtg = 6 at 
(a) 14 T and (b) 8 T. In the absence of equilibration and high electric field, resistance value for 
{i^tg,Vbg) = (6,-2) is different at different magnetic field. 


IX. EFFECT OF ELECTRIC FIELD ON MAGNETOTRANSPORT IN 
UNIPOLAR REGION 

In the unipolar region, the electric field is low compared to the bipolar region. So, to 
observe the modification of Landau levels by electric field in the unipolar region, we measured 
Rxx as a function of gate voltages from 2 T to 3.5 T. At a constant Vbg of 20 V, Rxx as a 
function of Vtg shows SdH oscillations as magnetic field is increased above 3 T as seen in 
Figure S12(a). 

We measured Rxx as a function of Vtg and Vbg at 3.5 T (Figure S12(b)) in the unipolar 
region. Line plot of Rxx vs Vbg at Ug = 1 V is shown by the green curve in Figure S12(c). 
At a constant B of 3.5 T, we observe Rxx to oscillate as function of Vbg for a given Vtg. 
The oscillations arise due to periodic decrease of back-scattering as the Fermi energy moves 
through different Landau levels in BG region. 

At Vtg = 2.8 V, the electric field is higher compared to Vtg = 1 V. So we took line slice of 
Rxx vs Vbg at Vtg = 2.8 V (orange curve in Figure S12(c)) and observed the fading of SdH 
oscillations. So with increasing electric field, the Landau levels come close to each other 
resulting in a decrease in the modulation in resistance. Thus the decrease in amplitude of 
oscillation is consistent with the modification of Landau levels with increasing electric field. 
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FIG. S12. Effect of electric field on resistance oscillations in unipolar region (a) Rxx as a function 
of Vtg at Vhg = 20 V at different magnetic field. SdH oscillations are observed when magnetic field 
is increased above 3 T. The plots at different magnetic fields are offset for clarity, (b) Rxx as a 
function of Vtg and Vbg at 3.5 T in the unipolar region, (c) Line plot of Rxx as a function of Vbg 
at Vtg = 1 V (green curve) and 2.8 V (orange curve). Amplitude of oscillations decreases with 
increasing Vtg and hence increasing electric held. 


X. MEASUREMENTS FROM ANOTHER DEVICE 

In this section we present measurements of another device at 14 T which showed similar 
magnetotransport behavior. This device has a low mobility of 2009 cm^/Vs, and the ratio of 
top-gate capacitance to back-gate capacitance is 11. The number of top-gates in this device 
is 35 and the period is 150 nm. 

Figure S13(a) shows the variation of Rxx as a function of Vtg and Vbg at 14 T. All the 
plateaus are not well resolved due to low mobility. Dissipationless transport is observed for 
{vtg,i^bg) = (2,2), (-2,-2) and (-6,-6). The maximum resistance in this device, observed in the 
bipolar region, is 58 kQ. 

Figure S 13(b) shows line slice of Rxx as a function of Vbg at Vtg of -2. The black dashed 
lines correspond to edge state equilibration for a single top-gate. The experimental curve (in 
green) coincides with the black dashed lines at {vtg-,i^bg) = (-2,-6). Thus this device having 
35 top-gates shows the plateau in Rxx corresponding to a single top-gate in the unipolar 


region. 
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FIG. S13. Charge transport measurement of another device in the presence of a magnetic field of 
14 T. (a) Rxx as a function of Vtg and Vhg. (b) Rxx as a function of Ubg for i^tg = ~2. Resistance 
plateaus at {vtgjVtg) = (-2,-6) and (-2,2) for a single top-gate are marked by black dashed lines. 
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